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Symbols 

A(z), B(z) Z-dependent amplitudes of two modes 

B static magnetic flux, gauss 

CART cloud and radiation testbed 

CCW counterclockwise 

CW clockwise 

D diameter of telescope 

£>diff minimu m diameter of total energy 

DIAL differential absorption lidar 

d distance between atoms 

d length of medium, cm 

E h E 2 transverse mode patterns 

FWHM full width, half maximum 

/ focal length of telescope 

K spatial frequency 

k constant between 1 and 1.5 for wavelengths in visible spectrum 

k (} propagation number 

L length 

LASE Lidar Atmospheric Sensing Experiment 

lidar light detection and ranging 

N(R) water vapor concentration between R l and R 2 

N b number of counts/sec from background solar light 

N d number of counts/sec from detector dark current 

N s number of counts/sec from water vapor DIAL signal 

n refractive index 

n e ff effective refractive index 
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/?0 original fiber refractive index 

R, r defined by equations (19) and (18), respectively 

R h R 2 ranges 

S/N signal-to-noise ratio 

STP standard temperature and pressure 

Soft lidar signal return amplitude at off-line wavelength 

S on lidar signal return amplitude at on-line wavelength 

T measurement tune 

V factor of proportionality (Verdet constant), min of arc/gauss/cm 

x, y, z coordinates 

|3 propagation constant 

P i propagation constant of forward progagating mode 

P 2 propagation constant of backward propagating mode 

T phase mismatch 

y line width, pm 

An amplitude of change in refractive index compared with n Q 

AX wavelength separation between laser on- and off-lines 

Act differential absorption cross section between on- and off-line wavelengths 

K coupling coefficient 

A spatial period of grating 

X wavelength 

Xq wavelength of incident light 

CT cross section 

<)),9 polar coordinates 
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1. Introduction 


1.1. Atmospheric Water Vapor 

Atmospheric water vapor exists mainly in the troposphere, i.e., the part of the atmosphere that is 
closest to the Earth, from ground to 10 km. The concentration decreases nearly exponentially with 
height, as shown hi figure 1 (ref. 1), and varies rather smoothly with longitude but decreases rapidly 
poleward. Even though water vapor only exists in small amounts, it plays a fundamental role hi many 
atmospheric processes such as the atmospheric energy budget, global water cycle, atmospheric chemis- 
try, and transport of pollution. 

Climate change is caused by complex hiteractions among many elements such as radiation, clouds, 
aerosols, precipitation, atmospheric circulation, and greenhouse gases. Water vapor is the main factor in 
cloud formation (which hnpacts precipitation) and in severe storm development. Water vapor is also the 
most important greenhouse gas and is responsible for mainta inin g the current clhnate in that it absorbs 
the longwave reradiation from the surface of the Earth. 

The surface temperature of the Earth depends on the absorbed solar radiation and the rate at which 
energy is reradiated from the surface of the Earth to space. The rate of reradiation in the atmosphere 
depends on the amount and the vertical distribution of clouds, aerosols, and greenhouse gases, espe- 
cially water vapor. It acts as a buffer for the rise and fall of temperature. When the temperature drops to 
the dew point, water vapor condenses, releasing heat, and the temperature rises again. On the other hand 
when the temperature rises, water vapor evaporates, taking up heat, and in this way cooling the 
atmosphere. 

Most climate models show that an initial warming of the climate results in a significant increase in 
water vapor that then amplifies the initial warming. Therefore an increase in other greenhouse gases, 



Figure 1. Water vapor concentration profiles for indicated latitudinal bands and seasons. (See ref. 1.) 



such as C0 2 , causes a wanning of the climate that causes an increase in atmospheric water vapor, which 
results hi further wanning when the water vapor absorbs reradiation from the Earth to the sky. The satu- 
ration vapor pressure increases exponentially with temperature; therefore, the abundance of water vapor 
can provide a very important increase hi atmospheric temperature. 


The latent heat of water vapor also channels energy into the atmosphere, and it is the main energy 
source for hurricane development. Understanding the distribution of water vapor around a hurricane 
may one day allow the prediction of hurricane intensity and direction. 


A better understanding of the horizontal and vertical distribution and transport of water vapor is 
essential for understanding the processes that regulate climate and to better predict future climate 
changes and their regional and global impacts. The use of lidar onboard NASA research aircraft has 
allowed the measurement of the atmospheric water vapor distribution during different campaigns 
around the world. In the future, lidar systems will be deployed on Earth-orbiting satellites to measure 
the global water vapor distribution (refs. 2, 3, and 4). 

1.2. Lidar Measurement of Water Vapor 

Lidar (light detection and ranging) is a method that has been used since the 1960’s to detect parti- 
cles or gases in the atmosphere. A laser is used to send out short pulses of light. The pulse travels with 

o 

the speed of light, 3x10 m/s. Some of the radiation is reflected by the atmospheric molecules and 
aerosols and this back-reflected radiation is collected by a telescope. The distance traveled by the light 
to a point in the atmosphere is equal to the speed of light times half the round-trip travel time from the 
laser transmitter to the receiver. The radiation is analyzed, and the density of the atmosphere is then 
determined. Where the density is higher, a stronger backscattered signal is received. The atmospheric 
density is higher closer to the Earth; therefore, the return signal decreases with increasing distance from 
the Earth. 

Lidar is an active remote sensing method, which means that light is sent out actively as a laser pulse 
instead of using light from the Sun. This method gives more freedom because the system is not depen- 
dent on sunlight. There are several different ways of using a lidar system. It is easy to mount it on a 
truck, an airplane, or even a satellite so that the measurements can be made at many different places. A 
lidar system hi orbit around the Earth has many advantages, such as high resolution and complete Earth 
coverage. 

Lidar uses a single wavelength that scatters off atmospheric molecules and aerosols but caimot dis- 
tinguish between different atmospheric species. There are two main kinds of scattering, Rayleigh and 
Mie. Rayleigh scattering, also called molecular scattering, is scattering of light by objects that are small 
in comparison with the laser wavelength. The scattered flux density is inversely proportional to the 
fourth power of the laser wavelength. Mie scattering, also called aerosol scattering, is scattering by par- 
ticles that are larger than the wavelength of the radiation. The scattered flux density is inversely propor- 
tional to X , where k is a constant between 1 and 1 .5 for wavelengths in the visible spectrum. 


To measure the density of an atmospheric species such as water vapor, a more advanced lidar tech- 
nique is used. This technique is a special development of lidar called differential absorption lidar 
(DIAL) (ref. 4). 
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1.3. DIAL Technique 


Differential absorption lidar can provide long-range high-resolution measurements of atmospheric 
H 2 0. DIAL is based on the fact that different atoms and molecules absorb different laser wavelengths. 
In DIAL, two different wavelengths are sent into the atmosphere, one that will be absorbed by the parti- 
cle of interest such as water vapor and another wavelength that does not experience water vapor absorp- 
tion. Both wavelengths are near each other so that they experience the same atmospheric molecular 
scattering. By comparing the decay of the backscattered light from the absorbed wavelength with the 
decay of the backscattered light from the nonabsorbed wavelength, the concentration of water vapor can 
be determined as a function of range. 

The simplified DIAL equation for determining the average water vapor concentration N(R) between 
the ranges R\ and R 2 is given by (ref. 5) 


N(R) = 


1 _ 

2 A g(R 0 


*l) 




( 1 ) 


where >V on and S 0 g- are the lidar signal return amplitudes at the on- and off-line wavelengths, respec- 
tively. The lidar wavelengths are selected to be on the peak of the water vapor (such as 946.0003 nm) 
absorption and in an unabsorbed region off the peak (such as 945.9 mn), Ac is the differential absorp- 
tion cross section between the on- and the off-line wavelengths, and (R 2 - R{) is the range cell over 
which the average concentration is measured. Equation (1) assumes that the wavelength separation 
between on and off laser lines is small ( AX < 0. 1 mn) . In this case, aerosol backscatter and atmospheric 
extinction errors are negligible, and the laser lines are selected so that there is negligible interference 
from other atmospheric gases. 

Figure 2 shows the concentration of water vapor in the troposphere measured with different meth- 
ods (private communication from Syed Ismail et al., Langley Research Center). The solid curve is a 
Lidar Atmospheric Sensing Experiment (LASE) measurement which uses the DIAL technique. As the 
figure shows the DIAL measurement corresponds very well to the other methods such as Raman lidar 
and balloon sondes. 

In the future, a DIAL laser system will be placed on a satellite in space to continuously measure the 
atmospheric water vapor. Such a system would make it possible to predict the weather a month in 
advance by measuring the global distribution of water vapor. The water vapor lines of interest for a 
space-based DIAL system are the 940-mn water vapor band lines as shown in figure 3. These lines are 
of interest because of their very high absorption cross section (refs. 4, 5, and 6). 


1.4. Importance of Narrowband Receiver Filters in DIAL System 

In an Earth-orbiting satellite optical detection system, some of the sunlight that is reflected by the 
surface of the Earth and the atmosphere will always enter the telescope as noise. By reducing the optical 
bandwidth of the receiver, the noise is also reduced. To make good measurements, the signal must be 
larger than the noise or the signal-to-noise ratio should be greater than 1 . 

The signal-to-noise ratio for a photon counting receiver system is given by (ref. 5) 


s _ 

N VW s + 2<JV i , + Af, ( ) 


( 2 ) 
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Figure 2. Concentration of atmospheric water vapor (grams of water vapor per kilograms of air) in the 
troposphere. 
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Figure 3. Water vapor absorption spectrum showing water vapor absolute lines A, corresponding absorption cross 
sections < 3 , and FWHM line widths Y- 
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To improve S/N, the signal return N s could be made larger by building a higher energy laser on the 
satellite, which is a very expensive option. If the measurement time T is increased, S/N increases but the 
horizontal resolution of the water vapor measurement decreases. To get good measurements, a compro- 
mise has to be made between S/N and the horizontal resolution. The background noise signal N b and the 
detector noise N d should be as small as possible; N d depends on the detector chosen and to make it small 
a very good low-noise detector is needed; N b should also be as small as possible. To do this, an optical 
filter is placed in front of the detector. The fdters used today are typically 300-pm FWHM but this could 
be reduced with fiber Bragg grating filters to 20 pm FWHM. Figure 4 shows a future space-based water 
vapor DIAL system. Figure 5 (private communication from Syed Ismail, Langley Research Center, and 
Susan Kooi, SAIC) shows a simulation of such a system orbiting at an altitude of 400 km using a 
500-mJ laser energy at 10 Hz (940 nm) with a 1-m-diameter receiver telescope having a detector quan- 
tum efficiency of 40 percent. The horizontal water vapor resolution is 100 km and the vertical resolution 
is 1 km. Using different cross sections at different heights is necessary because the water vapor 
concentration decreases a lot with increasing distance from Earth. The difference in the two sets of 
graphs is the filter bandwidth, either 300 or 20 pm FWHM. The daytime operation filter transmission is 
30 percent for the 20-pm filter and is 50 percent for the 300-pm filter. Cross sections are given in units 
of x 10 -24 cm 2 . When using a 20-pm filter instead of a 300-pm filter, the random error decreases by 
about a factor of 2 from altitudes of 5 to 10 km. 

1.5. Research Objectives 

The primary goal of this research effort is to demonstrate and characterize the performance of an 
ultra narrowband optical fiber Bragg grating filter at 940 nm, which corresponds to the water vapor lines 
of interest for space-based DIAL measurements. A tunable diode laser is used to determine the band- 
width and efficiency of the filter and then these measurements are compared with theoretical calcula- 
tions. Laser radiation enters a 30-cm-diameter Newtonian telescope, where an optical fiber is placed at 
the focal point. Once the radiation is in the fiber, it is transmitted to a reflective fiber Bragg grating 
which acts as an optical filter. The efficiencies of each leg of this filter concept are measured so that 
overall filter performance can be determined. Suggestions for improvements in performance are then 
stated. 

2. Theory of Fiber Bragg Gratings 

2.1. Bragg Gratings 

Fiber Bragg gratings are based on the development of diffraction gratings. A grating is made up of a 
periodic reflective surface or structure. When light falls onto the grating, wavelengths whose traveling 



Figure 4. Space-based water vapor DIAL system. 
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Figure 5. Modelling of error in water vapor atmospheric profile measured from space-based satellite for three 
different water vapor absorption cross sections with either 50-percent transmission, 300-pm FWHM filter or 
30-percent transmission, 20-pm FWHM filter. 



distance corresponds to an integer times the wavelength will constructively interfere making an interfer- 
ence maximum as shown hi figure 6. 

Radiation is incident on a crystal, where the distance between the atoms is d. The ray that is 
reflected from the lower layer of atoms must travel a distance 2 d sin 9 farther than the upper ray. For the 
reflected rays to reinforce each other, this extra path length must be an integer multiple of X, which is 
exactly what Bragg’s law tells us: 


nk = 2t/sin 9 (n = 1,2,3,...) 


( 3 ) 
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2.2. Fiber Bragg Gratings 

If a germanium-doped silica core fiber is exposed to ultraviolet radiation there will be a permanent 
change in the core refractive index. A periodic refractive index change (i.e., a grating) in a fiber is 
accomplished by exposing the germanium-doped silica core fiber to periodic ultraviolet laser light. 
Interfering laser beams are used to produce a mask that is then applied to the germanium-doped fiber 
resulting in refractive index changes where the beams interfere constructively and there is no change in 
the regions of destructive interference. The degree of index change depends on the amount of time the 
UV laser source is applied to the fiber. Thus very precise periodic index changes can be created in the 
fiber. 

When light propagates through a fiber where a grating has been induced and if the Bragg condition 
is satisfied, the forward propagating light will be strongly coupled to the grating, which results in light 
propagating in the backwards direction. To satisfy the Bragg condition, the difference in the propaga- 
tion constants of the two coupled modes must be equal to the spatial frequency K as 

Pj-(-p 2 ) = K = ^ (4) 

where P | and P 2 are the propagation constants of the forward propagating mode and the backward prop- 
agating mode, respectively, and A is the spatial period of the grating as shown in figure 7. The case of 
contradirectional coupling describes coupling when two modes with the same propagation constant are 
traveling in opposite directions. 

For contradirectional mode coupling, the total electric field can be described by 

E(x,z ) = A(z) E^{x) e + B(z) E 2 (x) e^ Z (5) 



Figure 7. Fiber grating and incoming and reflected beam modes. 
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By inserting equation (5) into the wave equation for modes in a fiber in the absence of any perturba- 
tion and solving these equations give the following two equations that describe how the amplitudes vary 
with z: 


II 

^ 143 

iTz 

kB e 

(6a) 

II 

, -iTz 
kA e 

(6b) 


where k is the coupling coefficient describing the coupling between the modes in the fiber and is there- 
fore dependent on the index, wavelength of operation, and the extent of the periodic perturbation; T is 
the phase mismatch given by the difference in propagation constants, (3j and P 2 , with the spatial fre- 
quency K subtracted; thus 


r = Pj - (-P 2 ) - AT (7) 

In our case, we require phase matching, implying T = 0 or 

Pj + P 2 = (8) 

Since the two modes are identical, the propagation constants are the same; thus 

P, - 02 - Y "*« (9) 

A '0 

where n e g- is the effective index of the mode, that is, the actual refractive index under these conditions, 
and 


Pi = P2 

k Q k Q 


where k 0 , the propagation number, is 


k 


0 


271 

A, 0 


( 10 ) 


( 11 ) 


and A 0 is the wavelength of the incident light. The spatial frequency K was defined in equation (4) as 



( 12 ) 


where A is the period of the grating. Equations (8) and (9) are used to get the following final expression 
of the period of the grating at the wavelength of operation (ref. 7): 
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A = 


(13) 


^■o 

2n eff 


2 . 3 . Derivation of Grating Reflectivity 

Equation (6b) is differentiated with respect to z and using equation (6a) we obtain 


d~B 2 d 
— = kB 

dz 


(14) 


When solving this equation and substituting the solution mto equation (6b), the following equations are 
obtained: 


B(z) = hj e K ~ + b 2 e ^ (15a) 

A(z ) = hj e K ~ -b 2 e K “ (15b) 

If we assume unit power incident in mode A on a periodic waveguide of length L, A(z= 0) = 1. There is 
no back-coupled wave beyond z = L; thus, B(z=L) = 0. Solving equations (15) with these boundary con- 
ditions, the following two expressions are obtained: 


-K L 

b = e 
1 2 cosh k L 

K L 

b 0 = 

2 cosh k L 

Thus, the equations describing the amplitudes of the modes are 

„ , . sinhK(z - L ) 

B(z) = ^ — — 

cosh kL 

, cosh k(z-L) 
cosh K L 


where B(z) is the reflected mode; therefore, we can define a reflection coefficient as 


B(z= 0) 
A(z= 0) 


-tanli K L 


(16a) 


(16b) 


(17a) 

(17b) 


(18) 


The energy reflection coefficient is then 


R = 



tanli” kL 


(19) 
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A medium of refractive index n with a periodic refractive index grating given by 


n(z) = n Q + An sin^-z 


( 20 ) 


has a coupling coefficient of 


it An 



( 21 ) 


Substituting equation (21) into (19) gives the final expression for the reflectivity of a fiber grating of 
length L as (ref. 7) 


R = tanh 


2 


n An L\ 


( 22 ) 


Figure 8 shows the grating reflectivity plotted from equation (22). As the length of the grating 
increases, the grating approaches 100 percent reflectivity for a given An. Also when the amplitude of the 
periodic refractive index change An increases, the reflectivity also increases. To get a grating with high 
reflectivity and narrow bandwidth, a compromise is necessary between the high An needed to get high 
reflectivity and the lower An needed to get a narrow bandwidth. The stars in the figure correspond to 
values measured in the laboratory and are discussed later. 
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Figure 9. Bragg grating bandwidth versus grating length at wavelength 946 mn. 


2.4. Bragg Grating Bandwidth Expression 

By solving the equation for the reflection coefficient (eq. (19)) for the first minimum at one side of 
the reflectivity peak and multiplying it by 2 to get the whole bandwidth, from minimum to minimum, 
we get the corresponding bandwidth of the fiber Bragg grating as (refs. 7 and 8) 


AX 


0 


^7<k 2 

nn efi'L 


7 7 1/2 

L + n ) 


(23) 


where Xg = X 0 is the wavelength that is reflected by the fiber grating. 

Figure 9 shows a plot of the Bragg grating bandwidth in pm from equation (23) versus the length of 
the grating. The wavelength used is 946 mn, which is also the wavelength that was used in the experi- 
ments. It shows that the bandwidth decreases when the length of the grating is increased. The bandwidth 
of the grating also decreases when An, the amount of the periodic refractive index change, decreases. 
The stars represent the values that have been measured in the laboratory and are discussed later. 


3. Experimental setup 


3.1. Typical Water Vapor DIAL Receiver Setup 

Figure 10 shows a general picture of how the fiber Bragg receiver would be used in a typical space 
DIAL application for the measurement of water vapor at 946 mn. The Newtonian telescope collects 
some of the light that has been backscattered by the particles in the atmosphere. This light is focused 
onto a single mode fiber that is connected to an optical circulator where light travels from port 1 to 
port 2. Attached to port 2 are fiber Bragg gratings which reflect light at a specific wavelength. The 
DIAL on- and off-line wavelengths at 946.0003 mn and 945.9 mn are reflected by the fiber Bragg grat- 
ings, enter the circulator again at port 2, and leave at port 3, where this light is detected by an avalanche 
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Figure 10. Fiber Bragg filter receiver concept for typical space-based DIAL application. 


photodiode. The signal from the detector is digitized and sent to a computer for storage and processing. 
The yokes are used to tune the fiber reflection wavelength. By stretching the fiber, the period of the 
grating is increased, and according to equation (13), a higher wavelength is reflected. 

3.2. Laboratory Experimental Setup 

To determine the optical bandwidth and efficiency of this receiver concept, the setup shown in fig- 
ure 1 1 was constructed. The laser used was an external-cavity tunable diode laser from New Focus Inc., 
Model 6238, which could be tuned from 930 to 966 nm, with a maximum power output of 30 mW and 
line width of 0.014 pm at 940 nm. The photodiodes are Thorlab DET210 high-speed photodetectors. 
The absorption cell (Infrared Analysis model number G-4- 10-PA) is a 20-m cell, which means that light 
is multiply reflected 20-m inside the cell. The optical circulator was made by Optics for Research, Inc., 
model number OC-3-FFT-946-X-Z. The fiber Bragg gratings made by Communication Research Centre 
Canada, have an on-line wavelength of 946.0 mn and an off-line one of 945.9 mn. The fiber used in the 
gratings is a single-mode Coming 1060 photonic fiber with a core diameter of 6 pm. The yokes are used 
to change the wavelength to be reflected by the grating. 

The output of the laser is split into two beams. One beam is detected by a photodiode that is con- 
nected to the oscilloscope and is used to monitor the output stability of the laser. The other beam is split 
once again, one beam goes into the water vapor absorption cell, while the other beam goes through the 
optical circulator unattenuated and enters the fiber at port 2 where fiber Bragg gratings are attached. The 
on- and off-line wavelengths are reflected by the fiber grating and they reenter the optical circulator and 
are then directed to port 3 where light is detected by a photodiode. The signal is connected to the oscil- 
loscope to monitor the two reflection peaks from the fiber grating. 

3.3. Optical Circulator Description 

The optical circulator works in two different ways depending on whether the light enters port 1 or 
port 2. Either unpolarized or polarized light from the receiver telescope enters fiber port 1 as shown in 
figure 12. The calcite crystal is a double refractive crystal, which means that it splits the beam into two 
beams that are orthogonally polarized. The rays that are polarized perpendicular to the optic axis (the 
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Figure 1 1. Bragg grating filter measurement setup. 



ordinary rays) pass through the crystal unaffected, but the rays polarized in the plane of incidence (the 
extraordinary rays) are shifted horizontally as shown. Light comes out of the crystal in two separate 
beams that are orthogonally polarized. The half-wave (see the appendix) plate rotates the polarization of 
the beams by 45° having a fast and a slow axis. The slow axis is inclined 22.5° to the horizontal and the 
fast axis is perpendicular to the slow axis. The half-wave plate rotates the component of the beam that is 
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polarized parallel to the slow axis by 180°; thus, the total polarization of each beam is rotated 45°. The 
Faraday rotator (see the appendix) rotates the polarization of the light back 45° so that when light leaves 
the Faraday rotator it is polarized the same way as before entering the half-wave plate. A Faraday rota- 
tor always rotates the polarization of light in the same absolute direction no matter from which direction 
light enters the rotator. The calcite crystal at port 2 lias the optic axis inclined at the same angle as the 
first crystal but in the opposite direction, which makes the two beams recombine into one beam. Light 
out of fiber port 2 is now coupled to a single-mode fiber which contains the fiber Bragg gratings. 

Light that is reflected by the fiber Bragg gratings enters again at port 2 as shown in figure 13, where 
the beam is split into two beams, o-beam and e-beam, which are perpendicularly polarized. The Faraday 
rotator rotates each beam 45°, rotating the polarization in the same absolute direction independent of the 
direction of propagation. The half-wave plate rotates the polarization 45° more. The Faraday rotator and 
the half-wave plate have rotated the polarization of the beams a total of 90°. The beam that was the 
o-beam (vertically polarized) in the first crystal is now the e-beam (polarized in the plane of incidence) 
and vice versa because of the total rotation of 90°. The e-beam is now shifted in the calcite crystal. Both 
the e-beam and the o-beam are reflected by the mirrors into the calcite crystal at port 3, which recom- 
bines the two beams into one beam. The unpolarized light that comes out of port 3 is detected by a 
photodiode detector. 

3.4. Telescope 

To find the focal spot size of the Newtonian telescope, the setup shown in figure 14 was con- 
structed. Laser light is coupled into a fiber that is placed at the focal point of a large parabolic mirror. 
Parallel light reflected by this mirror enters the receiver telescope consisting of a 30-cm-diameter, 
60-cm focal-length parabolic aluminum mirror and is focused at the focal point of the telescope. A 
charged coupled device (CCD) camera (Edmund Industrial Optics L55-702), pixels 537 (H) x 550 (V) 
and pixel size 7.15 (H) X 5.55 (V) pm, is placed at the focal point of this mirror to determine the focus 
spot size and is connected to an oscilloscope and also to a beam profiler which allows the beam spot size 
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to be viewed as measured. The focal spot size is limited by diffraction and theoretically the minimum 
diameter for 9 1 percent of total energy Z) diff is given by 

D„ - 4Axl (24) 

where /is the focal length of the telescope (60 cm), I) is the diameter of the telescope (30 cm), and X is 
the wavelength of the light (946 nm) which gives a minimum spot diameter of 8.3 pm. This focal spot 
size is close to the single-mode fiber to be used to couple the light into the optical circulator. 

4. Experimental Results 

4.1. Tuning of Fiber Bragg Grating 

When the fiber Bragg gratings were received from Communication Research Centre Canada, the 
bandwidth and wavelength of peak reflectivity were not known because the company did not have the 
means to measure these gratings. The first research step was to measure the wavelength of the peak 
reflectivity as well as the FWHM bandwidth. The setup shown in figure 10 was used to identify the 
reflection wavelengths of the fiber Bragg gratings. The tunable diode laser had its beam split into two 
beams, one coupled to the fiber Bragg grating using a Optics for Research model PAF-X-5-946 port and 
the other beam went to a water vapor absorption cell. The output of the cell had a silicon photodiode to 
measure the water vapor absorption lines. The diode laser was then tuned simultaneously through the 
fiber Bragg gratings and the water vapor absorption cell. A photodiode at the end of the Bragg grating 
fiber recorded the absorption and the bandwidth of the two Bragg gratings. 

As shown in figure 15, neither of the fiber Bragg gratings overlap water vapor absorption lines. The 
mask used to make these gratings was slightly off the designed wavelength of 946.0 nm for the on-line 
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Figure 15. Fiber Bragg grating reflection peaks and water vapor absorption lines. 



grating and 945.9 nm for the off-line grating. To adjust the gratings so that the on-line one coincided 
with the 946-mn water vapor absorption line and the other called the off-line did not coincide with any 
line, two aluminum yokes were constructed as shown in figure 16. Each fiber Bragg grating was 
epoxied to the ends of the aluminum yoke. The yokes consisted of an aluminum arc with a screw and an 
adjustment knob. By turning the adjustment knob, the screw pulled the arc apart and the fiber grating 
was stretched. 

By stretching the fiber Bragg grating, the period of the grating increases and higher wavelengths are 
reflected according to equation (13). The adjustment knobs were turned until the on-line grating coin- 
cided with the 946.0003-mn water vapor line and the off-line grating was tuned to 945.9 nm. The final 
result is shown in figure 17. 
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Figure 17. Final reflection peaks of fiber Bragg grating and associated water vapor absorption lines. 

4.2. Bragg Grating Reflectivity and Bandwidth Measurement 

The Bragg grating reflectivity and bandwidth were measured by using the setup of figure 1 1 . The 
tunable diode laser was scanned through the Bragg grating reflection peaks while the photodiode at the 
end of the fiber grating recorded two dips in the laser output as the laser scanned through the gratings. 
Knowing the photodiode voltage that corresponded to no grating reflectivity and knowing that if the 
grating were 100 percent reflective, then zero voltage could be measured by the peak reflectivity as 
shown in figure 18. As shown in the figure, the on-line grating has a FWHM of 26 pm and a reflectivity 
of 93 percent, and the off-line grating has a FWHM of 32 pm and a reflectivity of 9 1 percent. These val- 
ues should be compared with the theoretical value of equation (23) which gives a line width of 40 pm 
for a wavelength of 946 mn. The experimental value is a FWHM value, but the theoretical value is the 
distance between the first intensity minimum on both sides of the peak. This explains why the experi- 
mental values are slightly lower than the theoretical values. 

4.3. Telescope Coupling Measurements 

The setup shown in figure 14 was used to determine the spot size at the focal point of the telescope. 
The telescope has a diameter of 30 cm and a focal length of 60 cm. If parallel light is coming into 
the telescope, the minimum diffraction limited spot-size diameter at the focal point would be 8.3 pm. 
The typical minimum spot size was found to be 78 x 60 pm, which corresponds to an elliptical area of 
4680 pm 2 . This is a very large spot size compared with the area of the fiber to which it should be cou- 
pled. The diameter of the fiber is 9.3 pm, which corresponds to a mode field area of 67.9 pm 2 . Compar- 
ing the area of the fiber core with the spot size area results in an efficiency of 1 .4 percent. This should be 
compared with the measured coupling efficiency of the light at the telescope focal point compared with 
light emerging from the single-mode fiber. This efficiency was found to be 0.46 percent. 


17 




Wavelength 


Figure 18. On- and off-line fiber Bragg grating for wavelength of 940 nm. 


A multimode fiber was also tried in the experiment due to its larger diameter core. However it was 
found that a multimode fiber cannot be coupled to the circulator, the reflection peaks were smeared out; 
therefore, in this filter, a multimode fiber is not an option. 

4.4. Receiver Efficiency Measurements 

A setup like that of figure 1 1 was used to measure the efficiencies of this filter concept. The intensi- 
ties at all the different steps in the receiver were measured with a power meter. The efficiencies are an 
average of three measurements and are shown in figure 19. To maximize the transmission in the circula- 
tor, the alignment is crucial and quite complicated. At port 1 and port 2, x, y, z, 0, and (]) were adjusted 
carefully. It is necessary to iterate back and forth between the two ports to get the maximum output at 
port 2. The different efficiencies in figure 19 were measured step by step in the following way and they 
are summarized in table 1 . 

In measurement from path 0^1, light entered the telescope from the tunable diode laser. It was 
then reflected by the telescope into a single-mode fiber placed at the focal point of the telescope. The 
intensity before entering the telescope was measured compared with the intensity leaving the far end of 
the fiber. The efficiency was found to be 0.42 percent. 

Measurement from path 1— » 2 was made by comparing the measured intensity leaving the single- 
mode fiber before entering the circulator at port 1 with the measured intensity leaving the circulator at 
port 2, with the fiber port removed. The efficiency was 83 percent which is a very good value compared 
with the optimal value of 85 percent found by the manufacturer. 

Measurement from path 2 —> 4 was made by comparing the intensity leaving the far end of the fiber 
Bragg grating at a wavelength that was not reflected by the grating with the measured intensity leaving 
the circulator at port 2 with the fiber port removed. The efficiency was found to be 52 percent. 
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30-cm-diameter 
parabolic telescope 



laser Total efficiency 011 -line (0,1, 2, 4, 2,3) = 0.126 percent 
Total efficiency off-line (0,1, 2, 4,2,3) = 0.124 percent 

Figure 19. Receiver efficiency step by step. 


Table 1. Explanations for Figure 19 


Path 

Input measurement 

Output measurement 

Efficiency, 

percent 

Comments 

0-» 1 

Laser input at 946 mn to 
telescope 

Output from single-mode 
fiber before circulator 
port 1 

0.42 

Poor coupling of telescope 
spot size to single-mode 
fiber 

1 — > 2 

Output of fiber before 
entering fiber port 1 

Output of circulator at 
port 2, with port removed 

83 

Manufacturer gives opti- 
mal transmission of 
85 percent (without 
removing fiber port) 

2 — > 4 

Output of port 2 of circula- 
tor with port removed 

Output from end of fiber 
Bragg grating, laser tuned 
off grating 

52 

Poor coupling at fiber 
port 2 

2 -> 3 

Light sent into far end of 
fiber Bragg grating; input 
measured while fiber 
disconnected from port 2 

Output at port 3 

75 

Manufacturer gives opti- 
mal transmission of 
8 1 percent 


Measurement from path 2 — > 3 was made by using a separate fiber port to enter light at the far end 
of the fiber Bragg grating. The intensity leaving the fiber Bragg grating before entering at port 2 was 
compared with the intensity leaving the circulator at port 3 at a wavelength that was not reflected by the 
Bragg grating. The efficiency was found to be 75 percent, which is a very good value compared with the 
optimal value of 8 1 percent found by the manufacturer. 
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As can be seen, the greatest loss is the coupling between the telescope and the fiber. The efficiency 
of the filter without the telescope is 30. 1 percent for the on-line and 29.5 percent for the off-line. These 
numbers should be compared with the efficiency of a commercial 300-pm dielectric interference filter 
which is about 50 percent. Bear in mind that the bandwidth of fiber Bragg grating is much smaller. 

5. Discussion of Results 

When using optical filters in lidar receiver systems, the two major issues of primary concern are the 
optical transmission and bandwidth. Ideally the filter would have 100 percent transmission and a trans- 
mission bandwidth corresponding to the optical source bandwidth. The fiber Bragg grating described in 
this research effort approximates this ideal filter by having a reasonable optical transmission and very 
narrow bandwidth. It also has the ability to be easily tuned. 

To demonstrate and characterize the performance of the ultra narrow bandwidth fiber Bragg grating 
filter, an experimental program was undertaken to measure the bandwidth, grating reflectivity, and the 
overall efficiency of the filter. A reflectivity of 9 1 percent and a bandwidth of 32 pm were measured for 
the off-line grating (945.9 mn) and a reflectivity of 93 percent and a bandwidth of 26 pm for the on-line 
grating (946.0003 mn). The length of the gratings was 2.5 cm. The reflectivity can be made almost 
100 percent by making the grating somewhat longer, for example 3 cm; this would make the bandwidth 
slightly smaller. Another way of making the reflectivity higher is to choose a higher An, that is, a higher 
periodic change in the refractive index of the grating; this would be at the cost of a slightly broader 
bandwidth. 

The efficiency of the filter is about 30 percent after the telescope light has entered the single-mode 
fiber, which is good compared with commercial filters, which have an efficiency of approximately 
50 percent. Figure 19 and table 1 show that the transmission when passing through the circulator, port 1 
to port 2 and port 2 to port 3, is comparable with the transmission indicated by the circulator manufac- 
turer. The coupling from port 2 and into the fiber Bragg grating could be improved with better align- 
ment of the circulator; however, aligning the circulator is not trivial. The ports tend to get stuck and the 
alignment could be lost or could deteriorate with time or with changes in room temperature. There may 
be a better way of constructing the filter without using a circulator. The largest loss in the filter is when 
light is coupled from the telescope into the fiber. Theoretically the spot-size diameter should be limited 
by diffraction to 8.3 pm, which should couple well to a fiber with a core diameter of 9.3 pm. However a 
typical spot size measured was 78 X 60 pm, which corresponds to an area of 4700 pm 2 . If this area is 
compared with the core area of the fiber of 67.9 pm 2 , only 1 .4 percent of the light is coupled to the fiber. 
Experimentally an efficiency of 0.46 percent was measured. The difference could be explained by poor 
alignment of the fiber and the telescope. To get the minimum spot size at the focal point of the tele- 
scope, parallel light should enter the telescope. More research needs to be done on the coupling of the 
telescope light into the single-mode fiber. A larger fiber with a diameter of 68 pm was used, but it 
allowed multimode transmission and the circulator did not work correctly. The reflection peaks are 
smeared out, and one is not separated from the other. 

An important improvement in the signal-to-noise ratio can also be made with a fiber Bragg grating 
filter. Because the bandwidth is much smaller, less of the atmospheric background radiation and 
broadband light reflected by Earth is detected. An improvement of the signal-to-noise ratio by almost 
50 percent can be made as shown in figure 5. 

Figure 8 shows the reflectivity as a function of Bragg grating length for different values of refrac- 
tion index change An. With An of 5 X 10 -5 and a grating of 2.5 cm as was used in this experiment, the 
reflectivity should theoretically be 100 percent. Experimentally a reflectivity of 91 percent was 


20 



measured for the off-line grating and 93 percent for the on-line grating. The experimentally measured 
values are shown as stars in the figure. 

Figure 9 shows the bandwidth as a function of Bragg grating length for different values of refractive 
index change An. A 2.5-cm-long grating with a refractive index change An of 5 X 10 -5 , like the grating 
used in the experiment, should theoretically have a bandwidth of 40 pm. A bandwidth of 32 pm was 
experimentally measured for the off-line grating and a bandwidth of 26 pm for the on-line grating. The 
bandwidth calculated theoretically is the width between the first two m in ima on either side of the reflec- 
tion peak. The experimental bandwidth however is the FWHM, which explains why the experimental 
values are smaller than the theoretical one. 

Another advantage of this filter compared with commercial filters is that the on- and off-line wave- 
lengths can be tuned independently, which is not possible with the coimnercial filters. Rapid tuning 
could be accomplished by attaching the grating to a piezoelectric crystal. Commercial filters are too 
broad, covering both on- and off-line water vapor absorption lines. A fiber Bragg grating filter is 
compact and has low mass. It is therefore easy to transport and it would be possible to place it on an 
Earth-orbiting satellite. 

6. Concluding Remarks 

Fiber Bragg gratings have the potential for providing very narrow bandwidth optical filters for 
water vapor DIAL (differential absorption lidar) aircraft, space lidar receiver systems, and other filter 
applications. The filter bandwidth can be made as narrow as the water vapor absorption line, and sepa- 
rate filters can be used for on- and off-line DIAL atmospheric measurements. 

An ultra narrow bandwidth filter using fiber Bragg gratings was found to have a bandwidth of 
32 pm and a high reflectivity of about 9 1 percent for the off-line (945.9 nm) and a bandwidth of 26 pm 
and a reflectivity of 93 percent for the on-line (946.0003 nm). The measured values agree well with the- 
ory. Theoretically the reflectivity is 100 percent and the bandwidth 40 pm. However the bandwidth 
obtained theoretically is the bandwidth between the first two minima and the bandwidth measured 
experimentally is the FWHM (full width, half maximum). This difference explains why the experimen- 
tal bandwidth is actually smaller than the theoretical bandwidth. 

The on- and the off-line Bragg filters can be tuned separately, unlike coimnercial dielectric filters. 
The Bragg filters can be tuned easily by temperature or pressure. In the present research effort, the grat- 
ings were wavelength tuned by stretching them with an aluminum yoke. 

Overall filter efficiency of 30 percent was achieved after light was coupled from the telescope into 
the fiber as compared with 50 percent for coimnercial filters. However more research needs to be done 
on the coupling from the telescope into the fiber. 

An ultra narrowband fiber Bragg grating filter is a very promising way of measuring the water 
vapor content in the atmosphere from an Earth-orbiting DIAL satellite. With such a system, the under- 
standing of pollution and the prediction of weather and hurricanes would improve. 
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Appendix 


Description of Faraday Rotator and Half-Wave Plate 

Al. Faraday Rotator 

The Faraday effect is also called the magneto-optic effect. When linear polarized light propagates 
through a material medium such as glass and a strong external magnetic field is applied to the glass, the 
plane of polarization of the light will rotate an angle (3 (in minutes of arc) given by 

(3 = VBA (Al) 


where 

B static magnetic flux, gauss 

d length of medium, cm 

V factor of proportionality known as Verdet constant, min of arc/gauss/cm 

The Verdet constant for a particular medium varies with frequency and temperature. A positive Verdet 
constant corresponds to a material in which the Faraday effect causes a counterclockwise (CCW) rota- 
tion when looking at the light source, and the light travels parallel to the applied B field. If the light 
propagates antiparallel to the B field, it is rotated clockwise (CW) as shown in figure Al. A Faraday 
rotator is used when a rotation of the polarization of the beam is needed in the same absolute direction, 
regardless of the direction of propagation of the beam along its axis. This effect does not occur in the 
case of natural optical activity. By reflecting a beam back and forth through a Faraday rotator the effect 
can be amplified (ref. 9). 

The circulator hi our experiment used a Faraday rotator to undo the effect of the half-wave plate 
when traveling from port 1 to port 2. When traveling from port 2 to port 3, the Faraday rotator and the 
half-wave plate together changed the o-wave into the e-wave and vice versa by rotating the polarization 
of each beam 90° so that they were reflected by the mirrors and exit at port 3 . 



Figure Al. Faraday effect with V > 0. 
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A2. Half-Wave Plate 


A lialf-wave plate is made of a birefringent material, which means that the material has two differ- 
ent indices of refraction depending on the polarization of light: one for light polarized parallel to the 
optic axis of the material and one for light polarized perpendicular to the optic axis. The two polariza- 
tion states travel at different velocities through the crystal, with the effect that a relative phase shift is 
introduced. A half-wave plate introduces a relative phase difference of 180° between the two polariza- 
tion states. Light polarized parallel to the optic axis propagates with a higher speed than light polarized 
perpendicularly. Thus, the axis parallel to the optic axis is called the fast axis, and the axis perpendicular 
is called the slow axis. The slow axis in our experiment is inclined 22.5° to the horizontal and the fast 
axis is perpendicular to the slow axis. The half-wave plate rotates the component of the beam that is 
polarized parallel to the slow axis by 180°, thus rotating the total polarization of the beam 45° as 
explained in figure A2. Figure A2(a) shows the incoming polarization and its two components parallel 
to the fast and the slow axis. Figure A2(b) shows the outgoing light, after the component parallel to the 
slow axis has been rotated 180°. Figure A2(c) shows the incoming and the outgoing light, where the 
total polarization rotation is 45°. (See ref. 9.) 

The circulator in our experiment used a half-wave plate together with a Faraday rotator to rotate the 
polarization 90° when traveling from port 2 to port 3 . This rotation is done to change the o-beam into 
the e-beam and vice versa so that they are reflected by the mirrors and leave the circulator at port 3. 
When traveling from port 1 to port 2, the 45° rotation made by the half-wave plate is undone by the 
-45° rotation made by the Faraday rotator. 
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Optic axis 
or fast axis 



(a) Incoming polarization and two components of (b) Component parallel to slow axis rotated 180° 
incoming light. with new polarization of light shown. 


Optic axis 
or fast axis 



i 


(c) Incoming and outgoing light with total polarization rotation of 45°. 
Figure A2. Effect of half-wave plate. 
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